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ABSTRACT 

Current models of pulsar gamma-ray emission use the magnetic field of a rotating dipole in vacuum 
as a first approximation to the shape of plasma-filled pulsar magnetosphere. In this paper we revisit 
the question of gamma-ray light-curve formation in pulsars in order to ascertain the robustness of 
the "two-pole caustic" and "outer gap" models based on the vacuum magnetic field. We point out 
an inconsistency in the literature on the use of the relativistic aberration formula, where in several 
works the shape of the vacuum field was treated as known in the instantaneous corotating frame, 
rather than in the laboratory frame. With the corrected formula, we find that the peaks in the 
light curves predicted from the two-pole caustic model using the vacuum field are less sharp. The 
sharpness of the peaks in the outer gap model is less affected by this change, but the range of 
magnetic inclination angles and viewing geometries resulting in double-peaked light curves is reduced. 
In a realistic magnetosphere, the modification of field structure near the light cylinder due to plasma 
effects may change the shape of the polar cap and the location of the emission zones. We study the 
sensitivity of the light curves to different shapes of the polar cap for static and retarded vacuum dipole 
fields. In particular, we consider polar caps traced by the last open field lines and compare them to 
circular polar caps. We find that the two-pole caustic model is very sensitive to the shape of the polar 
cap, and a circular polar cap can lead to four peaks of emission. The outer-gap model is less affected 
by different polar cap shapes, but is subject to big uncertainties of applying the vacuum field near 
the light cylinder. We conclude that deviations from vacuum field can lead to large uncertainties in 
pulse shapes, and a more realistic force-free field should be applied to the study of pulsar high energy 
emission. 

Subject headings: MHD — pulsars: general — gamma-rays: theory — stars: magnetic fields 
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1. INTRODUCTION 

Pulsars are rotating neutron stars (NS) with very 
strong magnetic fields {B ~ 10^ — lO^^G). Nearly two 
thousand pulsars are known, mostly in the radio band. 
Compton Gamma Ray Observatory (CGRO) firmly de- 
tected seven gamm a-ray pulsars, and another three with 



less confidence [see Thompson ( 2004 ) for a review] . This 



number is rapidly increasmg with the start o f operations 



of Fe r mi Gamma-ray Space T elescope (e.g., |Abdo et al. 



2008 Abdo et al. 2009a|b|c |. Gamma-ray light curves 
dot 



of pulsars are typically double-peaked, generally out of 
phase with the radio, and have substantial off-peak emis- 
sion. Several theoretical models have been developed to 
explain the n ature of this emission, namely, the polar- 
ca p mod el ( Rudcrman & Sutherland '1975 Harding et 
al1|1978l paugherty & Harding.1982, .1996^ the sl ot-gap 



model (jjG, or two-pole caustic, TT^C tbr sh ort 
fc Scharleman n|[T979l pV rons"1983'; 'Muslimov 
2003 , 2004 ; Dyks fc Ruda k 2003; Dyks et al. 2004 
the outer-gap (OG) mode l (jCheng e t al. 1986a b) Ro 

Yadie 



Arons 
Harding 



and 



mani fc Yadigaroglu||1995 
In these models 



a 

tlie" 



20001 



garoglu.,1997j Cheng et 
particles are accelerated in 
regions, where strong electric fields are devel- 
to a deficit of charge. Gamma-ray emission 



gap" 
oped due 

is interpreted as the curvature and/or inverse Compton 
radiation from ultra-relativistic particles accelerated in 
these gaps. The models differ in the location of the gaps 
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in the magnetosphere. The polar cap model has narrow 
beam size and has difficulty in producing extended light 
curves, whereas SG/TPC and OG models can reasonably 
well reproduce double-peaked, extended light curves. Re- 
cent works have shown that OG model can also repro duce 
the high-energy spectrum for Crab and Vela pulsars ( Hi- 
roTani][2007| [Takata fc ChaKgl[2007i [Takata et aLl[20n7r 
2008| ). 

Special relativistic effects, such as the aberration of 
photon emission direction and photon travel time delay, 
are important for calculations of gamma-ray light curves 
when the emission zone extends far from the NS surface. 
In the TPC model, the emission zone is assumed to be 
along the last open field lines (LOFLs), extending from 
the polar cap to some cut-off radius. In the OG model, 
the emission zone is along the open field lines and extends 
from the null charge surface to the light cylinder^ (LC, 
i?LC — c/ft). In both models, the relativistic effects are 
essential to forming "caustics" in the sky map, which 
appear as sharp peaks in the light curve. The caustics 
arise when photons emitted from different regions of the 
magnetosphere happen to arrive to the observer at the 
same time. The presence and appearance of caustics is 
sensitive to both the relativistic effects and the geometry 
of the emission zone. 

All existing calculations of pulsar gamma-ray light 
curves assume that pulsar magnetospheric geometry can 
be represented by a vacuum dipole magnetic field. How- 

^ Recent version of OG model allows i nner boundary of em ission 
zone to shift toward the NS surface [e.g. [Takata et al.| | |2008[ l]. 
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ever, the pulsar masn etosphere is filled with plasraa ( Gol- 
dreich fc Julian|[l969 ), and plasma currents should mod- 
ity the held structure. The magnetosphere should then 
consist of the open and closed field line regions, sepa- 
rated by thin current sheets. In contrast, all field lines 
of the vacuum field, including those that travel beyond 
the LC, are formally closed, and no current sheets ex- 
ist. Numerical solutions of the structure of plasma-filled 
magnetosphere are now known in the limit of force-free 
(FF) MHD for axisym metric rota tors ( Contopoulos et al. 
1999| |Gruzinov^,2005| |Timokhin||2006, McKinney, 2006| ) , 
and, recently, tor three-dimensional oijlique rotators as 



well (Spitkovsky 2006 [Kalapotharakos &: Contopoulos 
2009 \. 'i'he J:''i'' field clearly demonstrates the current 



sheet structure, and its geometry differs substantially 
from the vacuum field near the LC [see |Spitkovsky^ (2008 j ) 
for a review] . Hence, the more realistic J^'i*' field geometry 
can lead to modification s of gamma-ray light curves [see 
Bai & Spitkovsky (2008| for preliminary results]. On the 
other hand, the vacuum field has been used over the years 
to obtain light curves that compare very favorably to the 
existing data. It can be argued that if the emission comes 
from the regions in the magnetosphere that are not too 
close to the light cylinder, the field geometry there may 
be well approximated by the vacuum field. This raises 
the question of how reliable and robust are the vacuum 
field light curves, and whether perturbations introduced 
by the presence of plasma in the magnetosphere would 
strongly affect the result. 

Although our ultimate goal is to study gamma-ray 
emission using the force-free field, in this paper we con- 
centrate on the modeling of light curves using vacuum 
magnetic field only. The results with the FF field will 



be pr esented in the companion paper ( Bai & Spitkovsky 
20101. We feel it is necessary to clarify a number of 



points before moving forward. As we try to reproduce 
the sky maps and light curves using vacuum field geom- 
etry, we find that there are ambiguities in the literature 
on the use of the aberration formula. Since the aberra- 
tion effect is crucial to the formation of caustics in any 
field geometry, in this paper we clarify the applicability 
of aberration formulas and compare their influences on 
the sky maps and light curves. In order to investigate 
the potential effects of plasma on the formation of light 
curves, we also study the sensitivity of the vacuum light 
curves to variations in the shape of the polar cap. Even if 
the field geometry in the bulk of the magnetosphere could 
be approximated by the vacuum field, it is the behavior 
of the field lines near the light cylinder that determines 
the shape of the polar cap and thus the location of the 
magnctosphcric emission zones. The plasma effects near 
the light cylinder can then undermine light curve mod- 
eling that uses the polar caps of the vacuum field. In 
this paper we show that the appearance of the sky map 
and light curves is indeed very sensitive to both the field 
geometry and the geometry of the emission zones, which 
suggests that revisiting theoretical models with the more 
realistic FF field is essential. 

This paper is structured as follows. We begin with 
the vacuum magnetic field formulas for pulsar magneto- 
sphere in section 2, and then discuss the effect of aber- 
ration in section 3. In section 4, we construct the shape 
of the polar caps. We present comparisons of sky maps 
and light curves between different aberration formulas as 



well as different polar cap shapes for the two-pole caustic 
and the outer-gap models in section 5. In section 6 we 
summarize our results. 

2. VACUUM MAGNETIC FIELD FORMULAS 

We will use the vacuum field as an approximation to 
the magnetic field of the plasma-filled magnetosphere, 
where E • B = and the magnctosphcric structure is 
stationary in the corotating frame (CF). Therefore, in 
the lab frame (LF), for any specified magnetic field B, 
there exists an electric field E such that: 



„ r2 X r 

E = X B . 



(1) 



where ft is the angular velocity of the NS, and r is the po- 
sition vector. Equation ([l| is of fundamental importance 
for discussing the aberration effect (§3), and is assumed 
throughout this paper. It does not apply, however, in the 
gaps where the high-energy emission is thought to origi- 
nate. We address this issue in Appendix B and show that 
the conclusions of this paper remain unchanged even in 
the presence of gaps. 

We will consider two commonly used formulas for the 
vacuum magnetic field, namely, the static dipole and the 
retarded dipole. For the static dipole, the field geometry 
is assumed to be the same as in a non-rotating dipole 
field, which is rigidly attached to the rotating pulsar. 
The field expression is 



B 



(2) 



where fl is the magnetic dipole moment vector and f is 
the radial unit vector. For a pulsar rotating along z axis 
with angular velocity and magnetic inclination angle 
a, the time evolution of magnetic moment vector is 

fl{t) — /i(sin a cos fltx + sin a sin flty + cos az) . (3) 

The static dipole is not a full solution of the field of the 
rotating dipole in vacuum. As such, it can be interpreted 
as known either in the corotating or the lab frame. As an 
approximation, we assume that the magnetic field of the 
static dipole is valid in the lab frame. The field geometry 
is sketched in Fig. [T] 

The full solution of the electromagnetic field of a ro- 
tating magne tic dipole is kn own as the retarded dipole 
formula (e.g., Jackson|[l975 ): 



B 



rr • 



,m , Jit) , m 



.(4) 

When r is small, the retarded dipole field configuration 
is almost the same as that of the static dipole. The devi- 
ation increases as the radius approaches the LC. One can 
find the e xpressions in Cartesian coordinates in |Cheng et| 
al. (20001. The retarded formula is valid in the LF, and 
not in the CF. 
Another vacuum formula that can be us ed to represent 



the magnetospher e is the Deutsch field (Deutsch 1955^ 
Michel fc Li|1999 ) , which includes corrections due to the 
finite size of the star R^. These corrections are of second 
order in Rn/RlCj and can be ignored for most pulsar 
parameters. Therefore, we will not consider the Deutsch 
field further^. 

^ In addition to the magnetic field, both the retarded dipole 



Modeling Gamma-Ray Pulsars using Vacuum Dipole Magnetic Field 



3 




Fig. 1. — Schematic illustration of magnetospheric geometry, a 
is the inclination angle of the magnetic axis, (^obs is the observer's 
viewing angle, and pQj = marks the surface of zero charge den- 
sity [pa J ^ -n-B/(27rc)]. 

3. RELATIVISTIC EFFECTS IN PULSAR RADIATION 

Pulsar gamma-ray radiation is believed to originate 
from curvature or inverse Compton emission from rel- 
ativistic particles (electrons and positrons). In the ex- 
tremely high electromagnetic field, these particles travel 
along magnetic field lines in a frame where electric field 
vanishes^ . Emission direction from these highly relativis- 
tic particles, therefore, coincides with their direction of 
motion. In order to collect the pulsar emission to obtain 
a light curve, we need to account for relativistic effects 
including the aberration of photon emission direction^ 
and the correction of photon travel time (time delay) . 

The time delay effect is straightforward to calculate. 
The difference of light travel time for photons emitted 
from different regions of the magnetosphere results in the 
delay in rotational phase at wh ich the photon is observed 
(Romani & Yadigarog hi][l995| 

(/-d = -r ■ v/Rlc , (5) 

where r is the position of the emission point, rj is the 
aberrated direction of photon emission (unit vector, in 
LF). The total phase of emission is = —(j)em + 4>d, where 
4>em is the phase of 77. In order to determine fj from any 
location in the magnetosphere, we need to consider the 
aberration effect. As the direction of particle motion de- 

and the Deutsch field solutions have the associated electric field. 
Retarded dipole solution satisfies E ■ B = 0, while the Deutsch 
field generally has E ■ B 7^ due to unipolar induction. Since 
we are using the vacuum magnetic field to approximate the FF 
field geometry, we discard these electric fields and use equation ([ij 
instead. 

^ In the force-free magnetosphere, E ■ B = [cf. eq.|l]l] and 
|E| < |B|. Therefore, at any point we can always find an instanta- 
neous frame in which the electric field locally vanishes. This is not 
assured, however, outside the LC for fields in vacuum. 

The aberration of light describes the change in the direction of 
photon propagation between different inertial frames. In the con- 
text of this paper, it refers to the correction of the photon emission 
direction relative to the direction of magnetic field in the LF. 



pends on the direction of the magnetic field, one has to 
be careful about the frame in which the field is known. 
The vacuum field formulas are in general valid in the 
LF as discussed in section 2. However, it was commonly 
implicitly assumed in the literature that these fields are 
valid in the instantaneous corotating frame, leading to 
discrepancies. For clarity, we will distinguish the follow- 
ing three frames: 

1. The lab frame (LF), which is the inertial observer's 
frame in which the pulsar is rotating around the z axis. 

2. The corotating frame (hereafter CF), which is a non- 
inertial frame that corotates with the pulsar and is the 
frame in which the field pattern is steady. It is related 
to the LF by the coordinate transformation 

X = X cos -V y sin , 

y' — —X sin fit + y cos fit , (6) 

z' = z, t' = t , 

where the prime denotes the coordinate in the CF. 

3. The instantaneous corotating frame (hereafter ICF), 
which is a local frame of an inertial observer instanta- 
neously moving at the corotation velocity. Therefore, it 
is defined only inside the light cylinder. 

We will discuss the aberration effect in these three 
frames separately. We pose the question as follows. The 
pulsar magnetic field configuration is known in the LF at 
t — and the field pattern corotates with the star. The 
electric field in the LF obeys equation ([T]). The photon 
emission direction at any position is along the magnetic 
field in a frame where the electric field vanishes. Our 
goal is to calculate the photon emission direction seen 
in the LF. We will do the calculation in three different 
frames and show that we reach the same result, as ex- 
pected. From now on, we add superscripts "C" and "/" 
to denote fields in CF and ICF respectively. Fields with 
no superscripts will always refer to fields in the LF. 

First, we do the LF calculation. Consider the motion 
of an emitting particle in the pulsar magnetosphere. The 
force-free condition requires E -I- /3o x B = 0, where f3o 
is the normalized velocity of the particle relative to c, 
or the direction of the emitted photon in the LF. This 
equation, when combined with equation ([I]), implies 

A) = /B + Aot , (7) 

where /^rot = x r/c is the normalized corotation veloc- 
ity, and / is a coefficient. For the emitting particle, we 
have \/3o\ — > 1, from which / can be determined. Solving 
equation ^ with |/3o| = 1 fixes the emission direction 
fj = We note that there are always two solutions 
associated with equation ([?]), one associated with parti- 
cles traveling along the magnetic field line in the CF, the 
other with particles traveling in the opposite direction^. 

Second, we consider the calculation in the CF. The 
relation of the electromagnetic field between the CF and 
the LF is ( ,Schiff.,1939, Gr0n.l984| ) 

= B , = E + /3rot X B . (8) 

^ At the LC, / = is an obvious solution where the particle 
corotates at the LC. This solution corresponds to the "backward 
moving" solution. However, we always pick the other solution, 
where the particle moves outward. 
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According to equation ([ij, we have E'-^ = 0. The photon 
emission in this frame is thus along B'-^ — B. We em- 
phasize here that CF and LF are related by a coordinate 
transformation rather than Lorentz transformation. In 
Appendix A, we show that the aberration formula calcu- 
lated from this frame has the same expression as equation 

^. . . . 

J"inally, consider aberration in the ICF. The magnetic 

field in LF and ICF are related by a Lorentz transforma- 
tion. Using equation ([!]) for electric field in the LF, we 
obtain 



where 7 = (1 



-Bt, B^ = Bp/7 , 



, 



(9) 



t)^^/^ is the Lorentz factor of corota- 
tion, Bt and Bp denote the toroidal and poloidal compo- 
nents of the magnetic field. Note that in ICF the poloidal 
magnetic field is smaller. Since E-^ = 0, the photon emis- 
sion direction is thus 



Vt 



I = ±BtlB'^ , 77^ = ±Bp/7i?^ , 



(10) 



where B'^ = ^ Bf + [1 — j3'^^^)B'^ is the total magnetic 

field strength in the ICF, and plus/minus sign corre- 
sponds to emission along/opposite to the magnetic field 
line. To get back to the LF, one should perform a Lorentz 
transformation to aber rate the direction of the photon 



77; the formula is ( jPyks fc Rudak||2003 1 

7(l + /3rof??0 



r/ + [7 + (7 - 



(11) 



Substituting equation (10 1 into equation (11) we find 
Bl ± B,B' 



Vt 



R2 

^0 



Vp 



±B'„ 



B2 



B 



p ' 



(12) 



where Bn — 



B^ is the total magnetic field 



strength in the LF. 

We note that although the ICF is defined only within 
the LC, equation (12) is valid anywhere. In fact, the 



ICF is only one special instantaneous frame in which E^ 
vanishes. We can construct other instantaneous frames 

is a gen- 
ie E X B 



where the electric field is zero, and equation ( 12 
eral result. For example, we can also choose t 
instantaneous drift frame (IDF). One can show that by 
transforming the LF field to IDF to calculate the aber- 
ration leads to the same result as equation ( [12^ . Simple 
algebra shows that equations (12) and ([7| are exactly 
equivalent as expected. Therefore, our proof is complete. 
In many of the previous studies (e.g., IRomani & Yadi 



:aroglu-1995';'Yadigaro glu"IT997l[CHe5"g et al.|:^0001|DyEs 

; Rudak 2003; Dyks et al.| |2004|, it was iniphcitly as- 
sumed that the vacuum field, given by any formula in 
section 2, is valid in the ICF instead of LF. The follow- 



ing equation was used instead of equation ( 10 ) 
r,l = ±Bt/Bo , 



Vp 



±Bp/Bo . (13) 
The corresponding photon emission direction is then 
PmtBo ± Bt B„ 



Vt 



Bo ± Plot Bt 



Vp 



j{Bo ± /3,otBt) 



(14) 



In the limit /3 — ?> 0, equations ( 12 ) and ( 14) agree. The 
differences between them are of the 0IdeI(D{/3f^^.) when 
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Fig. 2. — Polar cap shapes for the static (dashed) and retarded 
(thin sohd) dipole fields. The thick solid line is a reference circle 
corresponding to 9q = arcsin \/ Rn / RhC- 

Pmt *C 1- Since the difference between the two aberra- 
tion formulas scales as /3^Qt, then to first order in r/i?LC 
treating the vacuum di pole field as known i n ICF instead 
of LF, as was done in Dyks et al. (2004), is a reason- 



able approximation. However, as we shall see in §5, this 
introduces significant differences in the calculated light 
curves. In the other limit, /3i-ot — >■ 1, the two treatments 
of aberration result in different photon emission direc- 
tions. Interestingly, the backward solution of equation 



(12) agrees with the forward solution of (14), where the 
emission direction is perpendicular to the rotation axis. 

Since the vacuum field only approximates the FF field, 
one could choose to treat it as known in the LF, CF or 
ICF, as a matter of approximation. However, as shown 
in the next section, the last open field lines must be 
traced either in the LF or in the CF, but not in the 
ICF. This makes the treatment in a number of previous 
works not self-consistent. We will show in §5 that the 
two approaches can lead to substantial differences in the 
sky maps and thus the light curves. 

4. LAST OPEN FIELD LINES AND THE SHAPE OF THE 
POLAR CAP 

Determination of the last open field lines (LOFLs) that 
separate the open and closed regions in the magneto- 
sphere is needed to find the shape of the polar cap and 
the location of the emission zones. LOFLs are typically 
found by tracing field lines and checking whether they 
close inside or outside the LC. This tracing must be done 
either in LF or in CF, but not in the ICF. From equa- 
tion ([9]), the toroidal field in ICF coincides with LF or 
CF field, but the poloidal component is 7 times smaller. 
ICF field is thus ill-defined at the LC. As i? -> Rlc, 



B 



0, the field is purely toroidal, and no field line 



actually crosses the LC. Also, since ICF is not a global 
frame, local tracing of the field does not guarantee the 
tracing of a consistent set of field lines as seen from other 
frames due to aberration. Thus, the shape of the polar 
cap in ICF or in any other instantaneous frame is not 
well defined®. 

Another way to think about field lines is that these are trajec- 
tories of forcc-frcc particles in the magnetosphere. As seen in the 
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To find LOFLs, we use 3rd order Runge-Kutta integra- 
tion to t race the magnetic field lines in the LF. Bisection 
method (Dyks et al. 2004) is used to find the magnetic 



colatitude ^™ of the rim of the polar cap for every mag- 
netic azimuth 0„i . We show the resulting polar cap shape 
for both the static and retarded dipole fields in Fig. [2 
for magnetic inclination angle a = 60°. The polar cap 
shape of the retarded dipole as seen in Fig. [2] is com- 



monly found in the literature 
Cheng et al.|[2000| |Dyks et al.l 
works the retarded dipole field 



2004 



Yadigaroglu 1997 
Therefore, in these 
dipole field was traced in the LF, but 
for the computation of aberration the field was treated 
in the ICF, leading to an inconsistency. 

The vacuum field geometry appears similar to the FF 
field geometry near the NS surface, but the two are 



substantially different near the LC (Spitkovsky 2006). 
The polar cap shape, as determined by LOi"'Ls, is thus 
sensitive to the field structure near the LC. This may 
cause significant uncertainties in the geometry of the 
emission zones. As we show in the companion paper 
(Bai & Spitkovsky||2010 1, the polar cap shape in the FF 
field is, in fact, more circular and also larger than in 
the vacuum field. Thus, a circular polar cap may bet- 
ter characterize th e emission zone geomet ry when vac- 
Dyks & Rudak (12003) considered a 
a circular polar cap, given by 
In this paper, we also 



uum field is used 
static dipole field 



itF 



0rn = Oa = arcsin i/Rj^^^TRi]^. 
consider such a circular polar cap which is indicated as 
a reference circle in Fig. [2j We will refer to the polar 
cap found by tracing LOFLs as "traced polar cap," to 
distinguish it from the circular polar cap. 

In order to parameterize different field lines, we define 
open volume coordinates on the polar cap. The magnetic 
colatitude of the polar cap rim 6'™(0m) is generically a 
function of azimuth. For any point on the NS surface at 
the magnetic colatitude 9m and azimuth we define 
the open volume coordinate of this point to be (rov, 0m), 
where r^v = 9m/S'^[4'm)- Therefore, the rim of the 
polar cap and LOFLs correspond to Tqv = 1- AH open 
field lines have Tov < 1 while all closed lines have > 1 . 
We emphasize that open volume coordinates should be 
defined for magnetic fields in the LF or CF^. 

5. SKY MAPS AND LIGHT CURVES FROM VACUUM 
DIPOLE FIELD 

In this section we construct the sky maps and calculate 
the associated light curves for pulsar gamma-ray emission 
in vacuum fields. We consider both the two-pole caus- 
tic (TPC) and the outer-gap (OG) models for static and 
retarded dipole fields. In each case, we compare the re- 
sults using different aberration formulas [i.e., equations 
(12) and (14)]. We also compare the results using differ- 
ent polar cap shapes. 



CF, the electric field is zero, and the trajectory just traces the field 
line. In the LF, eq. ([7l determines the trajectory, including the 
rotation of the field paFtern and motion along the field line. Since 
B*-^ = B, these two viewpoints are consistent with each other. The 
trajectory in the ICF, besides being ill-defined at the LC, does not, 
in general, trace the same field line as in the LF or CF. Additional 
discussion of particle tra jectories can be found in Appendix B of 
|Bai fc Spitkovsky|poTo| l. 

' The polar cap shape of a retarded dipole has a "notch" (Fig. 
[2]|. A special treatment of open volume coordinates in this region 
was given in Dyks et al. (2004). We do not use this method, as the 
differences in the sky maps due to better resolution of the notch 
region are small. 



The TPC model is an extended version of the slot-gap 
(SG) model, where emission has until recently been as- 
sumed to come from a thin sheet centered on the LOFLs, 
i.e., = 1.* The emission zone extends from above the 
polar cap along the LOFLs up to a certain cut-off radius, 
where emissivity is assumed to drop to zero. The cut-off 
is described by r,„ax, the distance to the center of the 
NS, and i?max, the cylindrical radius to the rotational 
axis. These two parameters constrain the extent of the 
emission zone. For the OG model, the emission zone is 
assumed to come from a layer in the open field line region 
beyond the null charge surface (NCS) where Goldreich- 
Juhan charge density equals to zero^ [poj — -6/(2710)]. 
We adopt r^^ = 0.9 for the center of this layer. The emis- 
sion can extend to -Rmax ~ Rhc, and no extra cut-off is 
needed. In all the cases, emissivity is assumed to be 
constant along the field lines in the emission zone. For 
different field lines, emissivity is weighed by a Gaussian 



function centered at r^y, with width a — 0.025 (Dyks et 



al. 



2004|). 

^'he radiation from the emission zones, when projected 
to the observer as the NS rotates, produces the light 
curve. It is convenient to collect all photons from the 
emission zones as they fall on the sky and plot the inten- 
sity as a function of stellar phase (j) and the observer's 
viewing angle ^obs- The light curve is then obtained by a 
cut trough this sky map at the observer's viewing angle 

^obs- 

Mathematically, the sky map is a map from the emis- 
sion zone to the sky coordinates. At fixed Tov, the emis- 
sion zone is a two dimensional manifold, which can be 
parameterized by {(j)m,l)j where I is the length of the 
magnetic field line starting from the polar cap region, 
and d/dl is along the direction of the magnetic field in 
units oi Rlc- With constant emissivity along field lines, 
the Jacobian determinant of this transformation deter- 
mines the intensity on the sky map: 



,^obs) OC 



det 



9(0, Cobs) 



di(t)m,l) 



(15) 



When the inverse of the Jacobian matrix is singular, in- 
finite magnification is reached, and / — >■ (X). This means 
that light rays emitted from the neighborhood around 
the position {<j)m,l) arrive to the observer at the same 
time, which greatly strengthens the intensity. This is the 
analog of caustics in optics and strong gravitational lens- 
ing. In the case of the sky map, a "caustic" is often un- 
derstood as regions on the sky map with strong enhance- 
ment, where the determinant in ( 15 ) is much greater than 
1. The caustics on the sky map then correspond to peaks 
in the light curve. There is no guarantee of the existence 
of caustics on the sky map mathematically, and a general 
field structure or a variation in emission region geometry 
does not necessarily produce caustics. In the following 

* In the extended version of the TPC model, the emission zone 
is localized between rov = 1 and rnv = 1 — (5, where 5 describes 
the thickness of the gap ( [Venter et al.|[2b09[ l. We do not use this 
definition in order to directly compar e with earlier works (e.g., 
|Dyks fc Rudak|2003| jPyks et al.|2004^ . 

^ Recent developnient of OG niodel allows the inne r boundary of 
the e mission zone to be shifted inside the NCS (e.g., [Takata et al.| 
[2008 [ I . For our purpose, to compare different aberration tormulas 
and polar cap shapes, it suffices to consider emission beyond the 
NCS only. 
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subsections we will show examples of how the properties 
of caustics change as we use different field geometries, 
shapes of the polar caps and aberration formulas. 

5.1. Two-pole caustic model with static dipole field 

The T PC model with a static dipole field was consid- 
ered by Dyks & Rudak (2003). They used a circular 
polar cap, with 9m = arcsin yj Rn / RlCi rather than a 
traced polar cap. In the calculation of the sky map, they 
traced the field lines using the vacuum dipole formula Q . 
Therefore, the field was treated as in the LF, but equa- 
tion ( 14 ) was used for aberration, making the treatment 
not self-consistent. In this subsection we reproduce this 
result and compare with other sky maps obtained from 
self-consistent treatment using the static dipole field. 

In Fig. [Sj we show four sky maps and light curves us- 
ing the static dipole field with inclination angle a = 70°. 
In (a) and (b) we use the traced polar cap, while in (c) 
and ( d) w e use the circular polar cap. Aberration for- 
mula ( [l4| is used in panels (a) and (c), which is not self- 
consistent. Other tw o p anels have self-consistent a berra- 
tion [using form ula (12|]. The case considered by Dyks 
|fc: Rudak (2003) corresponds to Fig. |3|;. Representative 
light curves are plotted on the right ofeach panel in Fig. 
[3j the observer's viewing angle is chosen to be 80°. 

In all panels of Fig. [3j bright arc-like regions are caus- 
tics, which translate to sharp peaks in the light curve. As 
far as the bulk appearance of the sky maps is concerned. 
Figs. [3^ and[3]3 look similar, while Fig. [3j; andjSJi look 
similar. On the other hand, in each of the correspond- 
ing pairs, the details of the caustics differ. Therefore, 
the overall appearance of the sky map is very sensitive 
to the shape of the polar cap, while the choice of aber- 
ration formula is important for determining the detailed 
features of the light curves. To some extent, the incon- 
sistent treatment of magnetic field in the aberration for- 
mula in a number of earlier works may not be critical, 
but at the current level of very precise light curve mea- 
surement with the Fermi telescope, detailed differences 
in the light curves can be distinguishable. We now dis- 
cuss the sky maps and the corresponding light curves for 
the four cases. 

W e begin with Fig. [3};, wh ich is a reproduction of Fig. 



Dyks & Rudak (20031. In the sky map, we have 
two mam caustics near the polar caps, which are very 
prominent and sharp. Their ^obs extends from 50° to 
110° in the sky map. Right behind the main peak, a 
small hump in the sky map leads to a "subdominant" 
peak (labeled with "S"' in the light curves). As a result, 
two sharp peaks and a small, subdominant peak right 
behind the first peak is present in the light curve. This 
light cur ve was used to e xplain Vela's gamma-ray light 
curve by 
While 



Dyks fc Rudak] ( |2003| ). 
i^'ig. [5Jl looks similar to Fig. [3t 



a change in the 

aberration formula causes some quantitative differences. 
The "subdominant" peak becomes stronger, and the ex- 
tent of the main caustics becomes smaller. As a result, it 
is less probable to have a double-peaked light curve and 
more difficult to reproduce Vela's profile. 

Figs. |3|i,b, where LOFLs were traced to get polar cap 
shape, look quite different from Figs. J3b,d, which used a 
circular polar cap. Compared to Fig. [3p, the shape of the 
main caustic is more curved, and the original subdomi- 
nant peaks become strong and extended caustics. Also, 
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Fig. 3. — Sky maps (left) for the TPC model using static dipole 
field with inclination angle a = 70°. (a),(b) use traced polar cap 
while (c) (d) use circular polar cap. Aberration effect is t reat ed 
consistently in (b) and (d), while in (a) and (c), equa tion ||14|l is 
used f or aberration. Panel (c) is a reproduction of |Dyks &c Kudakj 
l |2003[ l. In all panels, lines are traced to rmax = O.yu. Dark circles 
indicate the polar cap. Color scale of each panel is independent, 
set by the minimum and maximum counts on the sky map. To the 
right of each panel we plot the light curves for observer's viewing 
angle ^obs = 80° . 

the emission zones resulting in two peaks now cover a 
smaller fraction of the sky map. Now it is almost impos- 
sible to reproduce Vela's light curve with this inclination 
angle. 

In sum, the double-peak feature for gamma-ray pulsars 
is hard to reproduce with the static dipole field, with the 
exception of case (c), where the inconsistent treatment 
of magnetic field in the aberration formula is used. This 
situation is true within a wide range of interesting incli- 
nation angles a (adjusting Tmax and i?max, or changing 
to smaller value does not help either). This exam- 
ple demonstrates how the sky maps and the resulting 
light curves are affected by the aberration formula and 
the shape of the polar caps. Meanwhile, the differences 
between the panels (b) and (d) in Fig. [3] characterize 
the uncertainties of sky maps and light curves calculated 
using the static dipole field. 
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Fig. 4. — Sky maps (left) for the TPC model using the retarded 
dipole field with inclination angle a = 60°. (a),(b) use traced 
polar cap while (c),(d) use circular polar cap. Aberration effect is 
treated consistently in (b) and (d), while in (a) and (c), equation 

fis used for aberration. Panel (a) is a reproduction of |Dyks et| 
|2004|| . Lines are traced to Tmax — 
1 and -Rmax = 0.75 m (aj, 
and (d); to rmax = 1.2 and i?max = 0.8 for (c). Dark circles 
indicate the polar cap. Color scales of (a) and (b) are the same 
to demonstrate the weakening of the caustic caused by aberration 
treatment. Similar for (c), (d). On the right of each panel shows 
the light curves at observer's viewing angle ^obs = 80°. 

5.2. Two-pole caustic model with retarded dipole field 

Next, we consider the t wo-pol e caustic model with a 
retarded dipole field (e.g., iDykTct al. 2 004) . We choose 
the inclination angle a = 60" and show the sky maps 
and light curves in Fig. [4j As in the last subsection, 
we also consider four cases. In all of them, the retarded 
dipole formula is considered valid in the LF, and is also 
traced in the LF. In (a) and (b), we use the traced polar 
cap, and the circular polar cap is used for the other two. 
In (a) and (c), the aberration formula assumes that the 
field is in the ICF (which is not self-consistent), while for 
the other two figures, the corrected aberration formula 
[eq. ([t])] is used. We have chosen ^obs = 80° for all light 
curves in Fig. [4j 

Fig. [4^ is a reproduct ion of the sky map and light 
curve of Dyks et al.| (2004). There are two strong caustics 




Fig. 5. — Sky maps (left) for the TPC model using the static 
dipole field with inclination angle a = 60° for comparison with 
Fig. |4] Circular polar cap is used. Aber rati on effect is treated 
consistently in (b), while in (a), equation | |14| l is used for aberra- 
tion. Lines are traced to rmax = 0.9. DarK circles indicate the 
polar cap. On the right of each panel shows the light curves at 
observer's viewing angle gobs = 80°. 



formed near the phase of two poles that extend over a 
wide range of observer's viewing angles. Just behind the 
first peak there is a small hump, caused by the overlap 
of emission from two poles. A small step-like drop in the 
mid dle corresponds to a weak discontinuity in polar cap 
rim ( |Dyks et aL||2004| . 

Fig! [3)3 is the corrected version of Fig. |4^, where 
aberration is treated self-consistently. Remarkably, the 
two strong caustics become blurred into weak and wide 
enhancements at the same phases. Since the caustics 
are weak, the brightest region in the sky map is mainly 
caused by the overlap of emission from both poles. In 
the light curve, we thus have two wide peaks, and a large 
fraction of radiation is from the off-peak phases. We have 
also explored the sky maps with other inclination angles, 
and find that the strong caustics appear in fewer com- 
binations of inclination angles and viewing geometries 
(see Appendix C]). The weak and wide caustic structure 
is generic for this TPC model using the retarded dipole 
field, and it has difficulty in producing sharp peaks in 
the light curve^°. 

Fig. |4j: and|4jl have circular polar caps. The appear- 
ance of the sky map is quite different from Fig. [4^ and 
|4]d. The caustics are more extended, and it is possible 
to get two peaks from the same pole (in total up to four 
peaks) . Meanwhile, we find the same trend: the inconsis- 
tent treatment of B field in the aberration leads to strong 
and narrow caustics; when aberration is corrected, the 
caustics become wide and weak. 

We see that the shape of the polar cap largely de- 
termines the overall appearance of the sky map, while 



If one changes r^^ to some smaller value (e.g., r^^ = 0.9), the 
caustics become less affected by the aberration formula. A full sky 
map from r°„ = 0. 9 using the retarded dipole field can be found in 
Fig. 12 of |Bai fc"Spitkovsky| | |2010| . The same trend can be found 
in the OG sky maps shown in §5.3. Although choosing smaller rJJ^ 
improves the light curves from the TPC m odel , it still suffers from 
large uncertainties as discussed further in 95.21 
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the details of the caustics are sensitive to the aberration 
formula. The geometry of the emission zone is mainly 
controlled by the shape of the polar cap. This example 
confirms that the shape of the polar cap (or the emission 
zone geometry) is one major source of uncertainty in the 
current modeling of gamma-ray pul sar light curves using 
TPC model (e.g., |Dyks etal][2004| . 

We can also directly compare the effect of different field 
geometries by using the same circular polar cap shape. 
Since we have previously chosen different inclination an- 
gles in Fig. [3] (static dipole, a = 70°) and Fig. [I] (re- 
tarded dipole, a = 60° ) i n orde r to reproduc e prev ious 



( |2004| , 



works by |Dyks fc Rudak| ( |2003 ); |Dyks et al.[ 
Figure [5] we plot the sky map for the static dipole field 
with a = 60° and a circular polar cap to facilitate the 
comparison with Fig. |4] We consider both treatments of 
aberration in Fig. [5] One can compare Fig. [4]:, |4ji with 
Fig. [5^, [5]d. The only difference between the two pairs 
of the sky maps is different field geometry. Since the re- 
tarded and static dipole fields are similar near the star, 
one might expect that the resulting sky maps would be 
similar as well. However, we see substantial differences in 
the appearance of the sky maps and the light curves, and 
the differences are already very prominent in regions that 
are not far from the star, as seen from the structure of 
the caustics. In particular, the main caustics form closer 
to the star and appear stronger in the case of the static 
dipole field than in the retarded dipole. Even though the 
caustics form at low altitude, they are very sensitive to 
the choice of field structure because caustics are a chance 
overlap of emission from different regions. This analysis 
shows that the appearance of the sky map is very sen- 
sitive to the field geometry itself, and the deviation of 
the vacuum field geometry from the more realistic FF 
field should be another major source of uncertainty in 
the modeling of gamma-ray pulsar light curves. 

5.3. Outer gap model with retarded dipole field 

The OG model with a r etarded dipole field was stud- 
ied by many authors (e.g., Romani & Yadigarog iu"l995! 
Yadigarogl ulflMTl |Cheng et al. 2000; Dyks et al. 2004, 
Takata fc C hang 2007| [Takata et al |2007[|2008l ). As be- 
fore, we consider four cases in J^'ig. 5: panels (a) and (b) 
with a traced polar cap, and c) and d) with a circular 
polar cap. Consistent treatment of aberration is done 
in (b) and (d) only. Case (a) was explor e d by Romani 
& Yadigaroglu ( 1995 ) ; Yadigarogiu] (|1997[) ; [Cheng et aT 



000|) ; iDyks et al.,(200l^ a nd case (b) was receii t lyex 
plore d by [ lakata fcT]IIang| ( [20071 ); [Takata et al.| ( |2007 

i'hese 



2008 1 , who have corrected t 



le aberration formula. 



studies were able to reproduce some of the gamma-ray 
pulsar light curves. Here we focus on the comparison be- 
tween the four cases to address the uncertainties of the 
OG model. 

Fig. |6] shows the sky maps for the four cases with incli- 
nation angle a = 65° . We further show the light curves at 
Cobs = 80° for all four cases to th e right of the sky maps. 
Our Fig. [6^ is a reproduct i on of [Romani fc Yadi garoglu 
■ r-^^ r^,^^ ^2000); Dyks et 



jl995); Y adigaroglu[ ( [1997| ; [Cheng et al 
[aLj ([2004| . I'he two peaks are caused by emission from 
the same pole. The offset from phase zero and the double 
peak profile are clearly present. The arrow indicates the 
effect of the "notch" in the polar cap, and the emission 
from both sides of the notch contributes to the caustic. 




Fig. 6. — Sky maps (left) of the outer gap model using the 
retarded dipole field at inclination angle a = 65°. {a),(b) use 
traced polar cap while (c),(d) use circular polar cap. Aberration 
effect is treated consistently in (b) and (d), while in (a) and (c), 
equation ( |14[ l is used for aberration. Panel (a) is the same as a 
number ot previous works. Lines are traced from the null surface 
to T-max = 1.6 and -Rmax = 1 for (a)-(d). The emission zone is 
centered at Tqv = 0.9. The arrows in (a) and (b) indicate the effect 
of the "notch" on the polar cap. On the right of each panel is the 
light curves at observer's viewing angle = 80° ■ 



In between the double peaks is the "bridge" emission 
with weaker intensity. In the plot with corrected aber- 
ration (Fig. |6]3), the regions on the sky map responsible 
for the first caustic become less extended. The emission 
from field lines on one side of the notch no longer forms 
caustics. Also, the regions responsible for the "bridge" 
become weaker. Fig. [6j: and [6 J are the sky maps for 
the OG model with circular polar caps. Although the 
sky maps look different from cases (a) and (b), the light 
curves look very similar to their counterparts [(c) resem- 
bles (a), and (d) resembles (b)]. 

We note that when we correct the aberration formula 
for the OG model, the sharpness of the two peaks is less 
affected than in the TPC model. This appears odd since 
the emission zone of OG is close to the LC, and one may 
expect larger difference between the two treatments of 
aberration. As one can see from Eq. (151, the intensity 
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of the caustics is a differential effect. It depends on the 
separation (on the sky map) between the projection of 
two neighboring points in the emission zone. Therefore, 
it is possible that caustic structures are similar for the 
two treatments of aberration in the OG model. 

Compared with the TPC model, the appearance of the 
sky maps and lights curves in the OG model is less sensi- 
tive to the shape of the polar cap. This is not surprising, 
because the open field lines not too close to the LOFL 
are not sensitive to small displacements in the polar cap. 
Therefore, the sky map does not change much when we 
replace the traced polar cap with a circular one. This 
might suggest that the OG model using the vacuum field 
is more robust than the SG model. However, the mag- 
netic field used by the OG model is needed close to the 
LC, which itself is still very uncertain. As we have dis- 
cussed in the last paragraph of §5.2, the uncertainties in 
the field geometry itself can also cause big differences in 
the sky maps. Therefore, the accuracy of the OG model 
using vacuum field is still in question. 

6. CONCLUSIONS 

All current models of pulsar gamma-ray light curves as- 
sume that the pulsar magnetosphere can be represented 
by the vacuum magnetic dipole field. The use of the 
vacuum field should be considered as an approximation 
to the field which includes the effects of plasma in the 
magnetosphere. In this paper we considered the applica- 
tion of vacuum field to theoretical models [the two-pole 
caustic (TPC) and the outer-gap (OG) models]. Our re- 
sults show that there are large uncertainties in using the 
vacuum field for predictions of gamma-ray emission. 

Calculations of high-energy pulsar light curves involve 
the construction of a map from the emission zone in the 
magnetosphere to the sky map. The appearance of the 
sky map is thus very sensitive to: 1) the geometry of 
the magnetic field; 2) the geometry of the emission zone. 
On the one hand, the vacuum field deviates from the FF 
field, which introduces one major uncertainty in the di- 
rection of emission. On the other hand, the shape of the 
polar cap is traced by the last open field lines (LOFLs), 
and sensitively depends on the field structure near the 
LC (Dyks & Harding 2004). This causes another large 
uncertainty because the geometry of the emission zone is 
largely determined by the shape of the polar cap. 

Relativistic effects including aberration and time de- 
lay are crucial to the formation of caustics in the sky 
map. We provide the treatment of the aberration effect 
in three different reference frames, namely, the lab frame 
(LF), the corotating frame (CF) and the instantaneous 
corotating frame (ICF), and show that they can be rec- 
onciled. To be consistent, however, the tracing of the 
field lines to find the polar cap must be done in the LF 
or the CF, but not in the ICF, because the latter is not 
a global frame. 

We compare the sky maps and light curves for the TPC 
model and the OG model using different aberration for- 



mulas. We find that the appearance of the caustics in 
the sky maps is sensitive to the treatment of B field in 
the aberration formula. For the TPC model using a re- 
tarded dipole field, we find that instead of having two 
strong caustics in the sky map, the corrected aberration 
formula weakens the caustics, leaving two wide and weak 
humps in the sky map. As a result, the conventional TPC 
model with retarded dipole field has difficulty in produc- 
ing sharp peaks in the light curve. For the OG model 
using a retarded dipole field and the corrected aberra- 
tion formula, the caustic responsible for the first peak 
in the light curve is less extended, and there is weaker 
"bridge" emission between the two peaks. Recent devel- 
opments in the OG model with the corrected aberration 
formula are still able to produce reasonabl e light curves 
( jTakata fc Chang||2007| [Takata et al.][2007| 

We study the uncertainties in the models of pulsar 
gamma-ray light curves by: 1) comparing the sky maps 
and light curves using different shapes of the polar cap, 
namely, the polar cap obtained by tracing LOFLs and 
the circular polar cap; 2) comparing the sky maps and 
light curves using the same polar cap shape (circular), 
but with different magnetic field configurations (static 
vs. retarded dipole). We choose circular polar cap as an- 
other possibility because t he polar cap of the FF magne- 
tosphere is more circular ( Bai &: Spitkovsky 2010) . Our 
results show that the overall appearance of the sky map 
is very sensitive to both factors. We find that for the 
TPC model, whose emission zone is centered on the last 
open field lines, the sky maps and light curves depend 
sensitively on the changes in the shape of the polar cap. 
Up to four peaks can be present when circular polar cap 
is used. This suggests that the repr o duction of the Vel a's 



light curve i nlDyks fc Rudak] ( [2003| ) ; |Dyks et al.| ( |2004| ) is 
not robust. For the OG model, whose emission zone is in 
the open field line region, we find that the sky maps and 
the light curves are not very sensitive to changes in the 
polar cap shapes. However, the vacuum field near the 
LC is unreliable, and the predictions from the outer gap 
model are still questionable. A detailed atlas of gamma- 
ray light curves for two models is presented in Appendix 
C. 

In all, we conclude that it is essential to revisit the ex- 
isting theoretical models using a more realistic magneto- 
spheric stru cture, i.e, the forc e-free field from numerical 

In the companion paper 



simulations ( Spitkovsky 



2006) 



Bai & Spitkovsky (12010), we will present the sky maps 



and light curves calculated using the force-free field. 
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APPENDIX 
A. ABERRATION FORMULA IN THE CF 
We can find the inverse transformation of equation ([6]) and write it in the differential form as 

fdx\ fcosM -sinm' -Oy\ /dx'\ 

dy' 

dz' ■ 
\dfj 




(Al) 

cdt' , we 
(A2) 



It suffices to consider t = t' = 0, in which case we also have x = x' , y = y' . Dividing each side by cdt 
obtain the photon velocity vector 

In this equation r/ is a unit vector that denotes the photon direction in the LF, while ff denotes the photon motion 
in the CF, but it is NOT a unit vect or. This is b ecause the metric in the CF is not the Minkowski metric, and has 
non-diagonal space-time components ( Schiff|[T939 1 . 

In the case of pulsar gamma-ray emission, ff is along the direction of B*^ = B. Therefore, equation (A2) can be 
rewritten as 

rf = /B-hO X r/c . (A3) 

This is exactly the same as equation ([t]). 

B. ABERRATION FORMULA IN THE PRESENCE OF GAPS 

In this appendix we relax the FF assumption to allow the presence of gaps in the magnetosphere. In the gaps, 
-Ey 7^ 0, and equation (|Tj) no longer holds. However, we show that under the following two assumptions, the aberration 
formula discussed in this paper still holds in the gaps: 

a) The magnetospheric structure is stationary in the CF; 

b) Particles move along magnetic field lines in the CF. 

Assumption 1) allows us to work in the CF using equation ([6|. Assumption 2) implies that in the CF, the perpendicular 
electric field must be zero (while E'^^ can be non-zero in the gaps). Transforming the field back to the LF by applying 
equation ([6]), we obtain 

r2 X r „ 

X B , 

(Bl) 



Ell 



E' 



CF 



Physically, assumption b) means the gap in the magnetosphere also corotates with the star, because the E x B drift 
velocity has a corotation component everywhere. 

Under the above two assumptions, even if particles are accelerated due to E\\ in the gaps, their direction of motion, 
and, hence, the direction of photons they emit, is always along the B field in the CF. Therefore, the treatment of 
aberration in this case is exactly the same as in the Appendix A, and the formula ( 12 ) is valid even in the presence of 
This analysis can be reproduced in the LF as well. 



gaps. 



The a bove discussions also apply if the magnetosphere rotates differentially with respect to the NS (e.g. 
2007a|b ), as long as one replaces f2 with the value of the angular velocity in the magnetosphere. 



Timokhin 



C. LIGHT CURVE ATLAS 



Knowing how the light curves a ppear for different view ing angles and magnetic obliquities can help place constraints 



on the pulsar geometry. Recently, [Watters et ah] ( |2009 ) compiled an atlas of TPC and OG light curves for a range of 
pulsar parameters using retarded vacuum field. 'I'he field was effectively assumed to be in the instantaneous corotating 
frame for the purposes of aberration calculation. Here we provide a similar atlas, emphasizing the differences brought 
by the self-consistent treatment of aberration. In Figs. [7] and |8] we show the sky maps and representative light curves 
for TPC and OG models, respectively. For every pulsar inclination angle (horizontal rows) we show the results for field 
taken in the lab frame (label LF, subpanels a, c, e, g) and in the instantaneous corotating frame (label ICF, subpanels 
b, d, f, h). As was shown above, only the LF treatment for the retarded vacuum field is self-consistent. The polar cap 
is obtained by tracing LOFLs in all cases. The light curve plots have identical vertical scales at the same obliquity for 
easy cross-comparison. 

The atlas displays the same general trends as seen before: self-consistent aberration reduces the strength of the 
caustics in the TPC model, and decreases the amplitude of the first peak for many observer angles in the OG model. 



Abdo, A. A., et al. 2008, Science, 322, 1218 
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Fig. 7. — Atlas of sky maps and light curves for the two-pole caustic (TPC) model for representative pulsar inclination angles (constant 
for each row) and for different observer viewing angles (arranged in columns). Panels a), c), e), g) show results from treating the retarded 
field in the lab frame, while b), d), f), h) treat the field in the instantaneous corotating frame. 
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Fig. 8. — Atlas of sky maps and light curves for the outer gap (OG) model for representative pulsar inclination angles (constant for each 
row) and for different observer viewing angles (arranged in columns). Panels a), c), e), g) show results from treating the retarded field in 
the lab frame, while b), d), f), h) treat the field in the instantaneous corotating frame. 
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